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A single magnetic monopole in pure SU(2) gauge theory is simulated on the lattice and its mass is computed 
in the full quantum theory. The results are relevant for our proposed realization of the dual superconductor 
hypothesis of confinement. 



1. INTRODUCTION 

In Ref. [1] we presented a realization of the 
dual superconductor hypothesis of confinement. 
The monopoles were of the 't Hooft-Polyakov 
(HP) form, in particular the Bogomol'nyi-Prasad- 
Sommerfield (BPS) solution. Assuming that 
the configurations considered dominate the string 
tension a we obtained the estimate y/a / A-g-^ = 
2.3. 

In our model we made two assumptions. The 
first one was that quantum fluctuations screen 
the (chromo-)electric field of the classical BPS 
monopole, turning it into what we call an HP-like 
monopole. The second assumption concerned the 
mass of the monopole. The mass of the classical 
BPS monopole is M = 4irfi/g 2 , where fi is the 
scale characterizing the solution. We conjectured 
that the mass of the monopole in the full quan- 
tum theory behaves as 
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Here g R is the running coupling in the i?-scheme 
defined in terms of the quark-antiquark potential, 
with A R /Aj^7 = 1.048 [2]. The function C was 
supposed to be slowly increasing as for the HP 
monopole, I < C 'S 2. 

The subject of the present study, described 



in detail in Ref. [3] (see Ref. [4] for preliminary 
work), is the assumption about the mass of the 
monopole. We put a single monopole background 
in a box and use lattice Monte Carlo methods 
to include quantum fluctuations. The mono- 
pole background is created by imposing appro- 
priate monopole boundary conditions. We use 
the simulation data to determine the function 
C(g 2 R (A R / n)) in order to discuss the validity of 
the mass formula (1). 

2. THE MAGNETIC MONOPOLE 

The classical magnetic monopole in euclidean 
pure SU(2) gauge theory is given by 
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as in the BPS limit of the HP monopole. The 
scale parameter /j, is arbitrary and can be re- 
garded as an inverse core size of the monopole. 
The configuration (2-5) is a static solution to the 
equations of motion, with mass M = Air/i/g 2 . 
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We want to study the monopole in a 4- volume 
consisting of a finite, non-periodic spatial box of 
size (2i?) 3 times a periodic time direction of ex- 
tent T. The idea is to induce the monopole by 
fixing the fields in the boundary at a value sug- 
gested by the asymptotic behaviour of the mono- 
pole field. We will assume /j,R > 1, so that effects 
exponential in /j,R are suppressed. 

The appropriate boundary conditions for the 
BPS monopole are given by its asymptotic be- 
haviour 

H(fir)/r ~ fi -!— , (6) 
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However, if our assumption is correct that the 
monopole becomes HP-like due to quantum fluc- 
tuations, 'HP-like' asymptotic behaviour, H/r ~ 
H, will be more appropriate at semi-classical val- 
ues of the coupling. 

We take HP-like boundary conditions with pa- 
rameter Ho f° r our dynamical simulations, i.e. 
K ~ 0, H/r ~ no- Although BPS-like and HP- 
like boundary conditions are not equivalent be- 
cause l/r(x) varies along the boundary of the cu- 
bic box, an analysis of classical monopole energies 
shows that our boundary conditions are compat- 
ible with both an HP monopole of scale \i = Ho 
and a BPS monopole of scale Heff = A*o + 1/ Ren- 

Another important result of the classical anal- 
ysis is that there is a symmetry implying that 
Ho boundary conditions are equivalent with 
monopoles characterized by scales h'o = A*o + 
2irn/T, for arbitrary integer n. This symme- 
try also allows monopoles of opposite electric 
charge, characterized by Ho < 0, to show up 
among the monopoles at positive values of Ho- 
All these monopole configurations are local min- 
ima of the action, but only the global minimum is 
important in the simulations. As a consequence, 
the accessible range of Ho values is restricted to 
< A*o < tt/T. We choose Ho = tt/7 1 to minimize 
the influence of nearby local minima. 

3. THE MONOPOLE MASS 

The monopole mass can be written as 
M = M m + M out , (8) 



with 



Z mon (j3\a,Ho) 



(9) 



^vac(/?) 

the contribution from inside the box, measured 
in the simulation, and M out a correction term for 
the outside region. Here 

Z mon {p-aHo) = j'DUexp[-S plaq (U;f3)] (10) 

is the partition function subject to monopole 
boundary conditions and Z vac (j3) is the analogous 
definition for vacuum boundary conditions, i.e. 
0. 

Eq. (9) can be written in a form accessible 
to Monte Carlo computation by differentiating it 
with respect to [3 and subsequently integrating 
again, 

MM = f^-AE(P). 

Jo p 

Here AE is the 'internal energy' 
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In order to compute M; n using Eq. (11) the 
integral is replaced by a sum. At each value of 
[3 in this summation two simulations are needed, 
to compute (S) with both monopole and vacuum 
boundary conditions. High statistics is required 
to compute the difference of the two large num- 
bers in Eq. (12) accurately. 

Fig. 1 shows the internal energy for a simula- 
tion at an 8 4 lattice, with uho = tt/8 = 0.39. 
In the low [3 region, AE « because the inte- 
rior of the box is decorrelated from the bound- 
ary. At intermediate [3 a monopole is induced 
in the box, and for large [3 the classical energies 
are approached. The interesting region extends 
from [3 = 2.3 onwards, where one expects scaling 
behaviour and the physical value of h increases 
from the order of the string tension scale to 
the high-momentum region. 

The AE data are integrated to obtain aM; n (/3). 
Subsequently, C 1Tl {g 2 R {KR / h)) is found as fol- 
lows. First Ar/h is calculated using Monte 
Carlo data for a^/a [5] and an input value for 
y/a I Kr (we used three different values [1,6,7]). 
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Figure 1. AE (in units of Air/ag 2 ) as a function 
of (3 = A/g 2 . Shown are the magnetic (O) and 
electric (A) components and the total (□). The 
horizontal lines denote the classical values. 



Next g^Aji/fi) is calculated using its two-loop 
/3-function and the corresponding C; n is extracted 
from aM m using Eq. (1). This procedure is car- 
ried out for both jj, = jj, and \i = corre- 
sponding to HP-like and BPS-like behaviour of 
the monopole, respectively. 

The result for C , including the correction for 
the exterior region, is shown in Fig. 2. For g\ w 
we know the monopole is BPS-like, so the HP 
analysis (upper three curves) is misleading there. 
At larger couplings, the monopole may or may 
not become HP-like. If it does, the C values will 
lie in the region indicated by the upper set of 
curves. This means that C increases from C = 1 
at weak coupling to C ~ 1.6 at g 2 R k, 8 or C ~ 2.0 
at g 2 R k, 10, depending on the value of yfa / Ar. 
This is in good agreement with our assumptions. 
If, however, the monopole remains BPS-like at 
large coupling, our first assumption does not ap- 
ply. Nevertheless, even in that scenario (lower set 
of curves) C remains almost constant, C ~ 1 up 
to g\ k, 6. 
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Figure 2. C as a function of g 2 R {AR/ fi) . The 
different sets of points are for y/a / Ar = 1.7 (O), 
2.0 (□) and 2.2 (A), for both jj, = jj, (upper set of 
curves) and /j, = (lower set). The horizontal 
lines denote the classical limit (g 2 = 0). The solid 
curves come from large-/3 fits of the mass data. 
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